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Abstract 

Using the DIANA data on the charge-exchange reaction K~^n —)> pK^ on a bound 
neutron, in which the s-channel formation of the pentaquark baryon 0'*'(1538) has 
been observed, we analyze the dependence of the background-subtracted 0''' —)• pK^ 
signal on the emission angle in the pK^ rest frame. In order to describe the 
observed cos 0™® distribution, invoking the interference between the nonresonant 
s-wave and the 0'*'-mediated higher-wave contributions to the amplitude of the 
charge-exchange reaction is required at a 2.8 cj level. The spin-parity assignment 
of 1/2“ for the 0'*' baryon is ruled out at a statistical level of 2.9 standard devia¬ 
tions. A physically-meaningful selection in cos 0™® based on the observed angular 
dependence of the 0'’“ —)• pK^ signal allows to boost the statistical significance of 
the signal up to 7.1 standard deviations. This is far in excess of previously reported 
signals and renders the 0"*“ existence more credible. 

PACS number(s): 13.75.Jz, 25.80.Nv 
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The exotic baryons with minimum quark configuration of (4g)g have been theo¬ 
retically discussed ever since the emergence of the quark model DEI. For such objects 
formed of light quarks, the lowest SU{3) representation was identified as the anti-decuplet 
that involves a single state with positive strangeness — the isosinglet baryon Q^{uudds). 
This pentaquark baryon can be uniquely identified by the KN decays {K^n and K^p) 
that are forbidden for the three-quark baryons. Since the “fall-apart” mechanism is not 
suppressed by any obvious selection rules, the decay width of the ©''■ baryon was phe¬ 
nomenologically assumed to be rather big (T ~ 100 MeV). The first rigorous predictions 
for the anti-decuplet of light pentaquark baryons were formulated in the landmark analysis 
[3] based on the chiral quark-soliton model. According to these theoretical predictions, 
the anti-decuplet baryons have spin-parity of and the mass of the isosinglet 0^ 

baryon should be close to 1530 MeV. The predicted decay width of the proved to be 
far below the earlier phenomenological estimates: T < 15 MeV. Subsequently, some the¬ 
orists using different assumptions came to a conclusion that the ©■*■ decay width should 
be well below this upper limit — on the order of 1 MeV or even less USE]. 

Narrow peaks near 1540 MeV in the effective-mass spectra of the systems nK^ 
and pK^ were initially observed in the reaction yu nK^K~ on the nucleus in 
the LEPS experiment [7], and in the charge-exchange reaction K^n —)■ pK^ on the Xe 
nucleus in the DIANA experiment [S]. Subsequently, both experiments confirmed their 
initial observations [91 uni dn [12]. Using the dynamics of s-channel formation of the 
©■•■ in the charge-exchange reaction K^n —)■ pK^, DIANA was able to directly probe 
the ©■*■ decay width; T = 0.34 ± 0.10 MeV assuming J = 1/2. Other searches for the 
©■•■ baryon in different reactions and experimental conditions yielded both positive and 
negative results, see the review papers [131 El E]. A number of experimental groups 
have reneged on their initial positive evidence, that anyway was statistically insignificant 
and may have resulted from wishful thinking and the so-called “bandwagon effect”. Of 
the many null results, only a few that have been formulated in terms of the ©■*■ intrinsic 
width should be treated as physically meaningful. The best (albeit mo del-dependent) 
null result has been reported by the E19 experiment at J-PARC, where the ©■*■ signal was 
searched for in the K~ missing mass in the hadronic reaction n~p K~X [Tn|. The E19 
upper limit on the O’*" decay width, P < 0.36 MeV assuming the ©■*■ spin-parity of 1/2''', 
is narrowly consistent with the DIANA measurement. On the other hand, a group from 
the CLAS collaboration has recently re-analyzed their data for the reaction yp K^K^p 
on hydrogen, invoking the interference between (f)p and ©’•'A'*’ in the final state pK^Kg 
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HT]. A narrow statistically-significant peak near 1540 MeV, tentatively interpreted as the 
signal, has been observed in the missing-mass spectrum. This observation does 
not contradict the null result earlier reported by CLAS for the same data sample [18] . 

In this paper, we continue the investigation of formation in the charge-exchange 
reaction K^n —>■ pK^ on a bound neutron using the data of the DIANA experiment. In 
particular, we probe the angular distribution of decay products in the G"'' rest frame. 



Figure 1: The cosine of the emission angle in the pK^ rest frame, 0^*^, for all mea¬ 
sured events (a) and upon applying the selections pt < 300 MeV and Pl > 0 (b). The 
corresponding distributions of simulated events are shown by crosses, and the contribu¬ 
tions of rescattering-free events — by dots. 

The DIANA bubble chamber hlled with liquid Xenon was exposed to a sepa¬ 
rated beam of monochromatic mesons from the 10-GeV proton synchrotron at ITEP, 
Moscow. In the hducial volume of the bubble chamber, momentum varies from ~ 730 
MeV for entering kaons to zero for those that range out through ionization. Throughout 
this momentum interval, all collisions and decays of incident mesons are efficiently 
detected. The momentum at interaction point is determined from the spatial dis¬ 
tance between the detected vertex and the mean position of the vertices due to decays 
of stopping mesons. Charged secondaries (electrons, pions, kaons, and protons) are 
identified by ionization and by decays at rest for kaons, and momentum-analyzed by their 
range in Xenon. The detection efficiency for y-quanta with p^ > 25 MeV is close to 
100%. Secondary mesons are identified by the detectable decays —)■ tt+tt” and 

—)■ 7r°7r°, and momentum-analyzed using the kinematic reconstruction. (In this anal¬ 
ysis, only the former decay is used.) Further details on the experimental procedure may 
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be found in [12] and references therein. The candidate events for the charge-exchange 
reaction K^n —>■ K^p with no intranuclear rescatterings are selected as hnal states with a 
single proton and a Kg —)■ 7r’''7r“ decay. The instrumental thresholds for the momenta of 
secondary particles are px > 155 MeV and Pp > 165 MeV. The experimental resolution 
is near 3.5 MeV for the pK^ effective mass. 

Plotted in Fig. [Tj^a) for all measured events is the cosine of the emission angle 
in the pK^ rest frame with respect to the pK^ direction of motion, cos0™®. Also shown 
is the cos0^® distribution of all simulated pK^ events (crosses) and of those in which 
the proton and the suffered no intranuclear rescatterings (dots). The former has been 
normalized to the number of all measured pK^ events. The effect of the selections in the 
transverse and longitudinal momenta of the pK^ system, pt < 300 MeV and > 0, 
is shown in Fig. [I](b). These are seen to reject the rescattered events rather than the 
unrescattered ones. The simulation procedure has been described in [T2] . 
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Figure 2: The mean value of | cos0™®| as a function of the pK^ effective mass for all 
measured events (a) and for those in the region 445 < Pbeam < 535 MeV (b). The effect 
of the selections pt < 300 MeV and > 0 is shown in (c) and (d). 
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The mean value of | cos0™’^| is plotted in Fig. |2]as a function of the pK^ effective 
mass. The enhancement observed at m{pK^) ~ 1540 MeV is emphasized by the selection 
in the K~^ momentum at interaction point, 445 < pbeam < 535 MeV, that reflects the 
dynamics of s-channel formation of the 0+ baryon in the reaction K~^n —?• pK^ on a 
bound neutron [12]. It is further emphasized by the selections px < 300 MeV and Pl > 0 
aimed at rejecting the rescattered events. That the anomaly in | cos0™®| occurs in the 
mass region of the observed 0+ peak [12] suggests that it is rooted in an “anomalous” 
angular distribution of 0"^ decays that may show a quadratic term in cos 0^®. Therefore, 
it is interesting to compare the cos 0^® distribution for the mass region of the peak with 
that for the sideband areas of m{pK^). Except for the scatter plots discussed in the last 
paragraph, the selection 445 < Pbeam < 535 MeV is implicitly assumed throughout. 




Figure 3: The pK^ effective mass prior to (a) and upon (b) applying the selections 
Pt < 300 MeV and > 0. Either mass spectrum is htted to the simulated nonresonant 
background with variable normalization plus a Gaussian with variable position, width, 
and magnitude. The null hts to the background form alone are shown by dashed lines. 


The distribution of the pK^ effective mass is shown in Fig.[3]^a), and upon applying 
the selections px < 300 MeV and p^ > 0 — in Fig. El^b). Either mass spectrum is 
then htted to the simulated nonresonant background with variable normalization plus a 
Gaussian with variable position, width, and magnitude. The width of the observed 0^ 
peak near 1538 MeV is consistent with the experimental resolution of cXm — 3.5 MeV. 
In agreement with the htted width of the 0^ signal, the peak area of the pK^ ehective 
mass is selected as 1530 < m{pK^) < 1546 MeV, and the sideband areas — as 1514 < 
m{pK^) < 1530 MeV and 1546 < m{pK^) < 1562 MeV. The cos 0™® distributions of 
events in the peak and sideband areas are shown in Figs. Hl^a) andlU^c), respectively. The 
ehect of the selections px < 300 MeV and p^ > 0 is shown in the corresponding right- 
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Figure 4: The cos0™® distributions for the G’*' mass region of 1530 < m{pK^) < 1546 
MeV (a) and for the sideband regions of 1514 < m{pK^) < 1530 MeV and 1546 < 
m{pK^) < 1562 MeV (c). The effects of the selections pr < 300 MeV and > 0 are 
shown in (b) and (d). The simulated distributions are depicted by dots. 


hand panels. For the sideband areas, the simulated distribution (dots) is normalized to 
the observed one by the number of events. Then, the same scaling factor is applied to 
the simulated cos0™® distribution for the peak area. (As a result, there the simulated 
COS0™® spectrum runs lower than the observed one.) The simulation that assumes a 
pure s-wave for the nonresonant reaction K^n —)■ pK^ |T9] agrees with the data for the 
sideband areas, but not for the peak area. 

In order to obtain the “pure” cos 0™® spectrum for the decay 0"^ —)■ pK^, we 
subtract the (halved) cos 0^® distribution for the sidebands from that for the G"*" peak 
region. Under the selections pt < 300 MeV and Pl > 0 that reject nearly a half of the 
rescattered events, the sideband-subtracted cosG^^® spectrum is shown in Fig. [5l It is 
then compared with Monte-Carlo predictions for the G"*" formation and decay assuming the 
angular distribution in the forms dlF/dcos 0^® ~ const, and dW/d cos ~ cos^ G^p® 
as shown in Figs. [5](a) and 151(b), where the simulated distributions have been scaled to 
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Figure 5: Under the selections pt < 300 MeV and pl > 0, the sideband-subtracted 
cos 0^® distribution for the 0+ mass region compared with Monte-Carlo predictions for 
the 0+ decay angular distribution in the assumed forms dW/d cos ~ const, (a), 
dW/dcosQ^^^ ~ cos^ 0™® (b), and dlU/dcos0™® ~ (cos0™® + + b with a = -I-0.2 

and 6 = 0 (c). The simulated distributions have been scaled to the data by area. 


the data by area. These forms are the two allowed components of the cos 0^® spectrum 
for arbitrary spin-parity of the decaying 0+ baryon. For both hypotheses, the values 
of x^/ndf (19.1/9 and 21.6/9) are unacceptably high. We have also verihed numerically 
that any linear combination of the above forms of the cos 0^® distribution for the decay 
0’'' —)■ pK^ leads to values in excess of 19 when compared with the observed sideband- 
subtracted cos 0^® spectrum. 

In order to verify the conclusions reached with the analysis of our low-statistics 
data, we also use the likelihood criterion. The likelihood function is constructed as a sum 
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—2 In L = 2 ^ [—rij + Ui + hi + rii In 


rii 


2 = 1 


Ui + hi 


where n* and z/j are the bin contents of the experimental and simulated cos 0^® distribu¬ 
tions for the 0’'' central mass region shown in Fig. 111(b), and hi are those of the simulated 
cos0^® distribution for the given hypothesis of 0^ decay. The latter distribution has 
been scaled by the number of events to the difference between the former two, so that the 
htting function as a whole is normalized to the observed cos0^® distribution by area. 
For the hypotheses dW/d cos ~ const, and dW/d cos ~ cos^ 0^h we obtain 
—21nL = 24.2 and 23.1, respectively. These values of —21nL for ndf = 9 correspond 
to the p-values near 0.005 and .007. As with the x^ analysis above, we hnd that assum¬ 
ing the 0’'' decay angular distribution in the form of an arbitrary linear combination of 
the former two fails to tangibly reduce the value of —21nL. This simple analysis sup- 
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ports the above conclusions based on the criterion. From the derived p-values we 
may conclude that the uniform angular distribution dW /~ const, is inconsis¬ 
tent with the data at a statistical level of 2.9a, and the more general symmetric form 
dhF/dcos©™® ~ a -f 6cos^ 0™® — at a slightly lower level of 2.8a. 
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Figure 6: The pK^ effective-mass spectrum upon rejecting the events that fall within 
the COS0™® intervals of —0.6 < cos05P® < +0.2 (a), —0.7 < cos©™*^ < +0.3 (c), and 
—0.75 < COS0™® < +0.35 (e). The effect of the selections pt < 300 MeV and > 0 
is shown in the right-hand panels (b), (d), and (f). Each mass spectrum is htted to 
the simulated nonresonant background with variable normalization plus a Gaussian with 
variable position, width, and magnitude. The null hts to the background form alone are 
shown by dashed lines. 















































The disagreement with any viable form of the 0^ decay angular distribution is 
rooted in the marked forward-backward asymmetry of the observed cos spectrum. 
The asymmetry in the form of a linear term in cos0^®, required at a level of 2.8a, may 
arise only from the interference [20] between the nonresonant s-wave and the ©■''-mediated 
higher-wave contributions to the amplitude of the charge-exchange reaction K~^n —)■ pK^. 
However, the interference should not affect the cos 0^® distribution for the signal as 
soon as the ©''■ spin-parity is 1/2“ implying an s-wave decay and a uniform cos0^® 
distribution. Therefore, the 1/2“ assignment is ruled out by the data at a level of 2.9 
standard deviations. 

The shape of the background-subtracted ©"^ signal in Fig. O suggests an angular 
dependence of an asymmetric quadratic form dW/d cos 0™® ~ (cos 0™® + a)‘^ + b, where 
6 ~ 0 and the offset parameter a is positive at a 2.8a level. Tentatively substituting 
6 = 0 and a = +0.2, we obtain y^/ndf = 4.3/9 with the approach as shown in 
Fig. El(c), and —21nL = 6.1 for ndf = 9 with the likelihood criterion which corresponds 
to a p-value near 0.8. (The assignments a = -1-0.1 and a = -1-0.3 also yield acceptable 
values of —21nL.) The detailed interpretation of this form of the background- 

subtracted cos 0^® distribution requires a theoretical analysis of interference effects in 
terms of helicity amplitudes, and therefore is beyond the scope of this paper (apart from 
excluding the ©''■ spin-parity assignment of l/2“). Instead, our major objective is to 
formulate a physically-reasonable data-driven selection that may render the 0^ signal 
more significant. This is an important task since the signals reported thus far iiiniini 
are but slightly in excess of 5cr, whereas the “credibility threshold” for proving the 0^ 
existence has been estimated as 7a [21] given the controversial experimental situation. 

If formation of the 0^ baryon indeed follows the distribution dW/d cos ~ 
(cos0™® -|- a)^ -|- 6 as argued above, rejecting the events with cos0™® values near the 
minimum of the parabola at —a should enhance the signal-to-background ratio in the 
pK^ effective-mass spectrum and the statistical signihcance of the ©''■ peak. Shown in 
Fig. El are the effects of cutting away the cos0^® intervals centered on cos05f'® = - 
0.2 : -0.6 < cos0^® < +0.2, -0.7 < cos0g^® < +0.3, and -0.75 < cos0^® < 
+0.35. Despite the uncertainty of the angular-distribution parameters, these selections 
are physically meaningful rather than arbitrary. Each mass spectrum is again fitted to 
the simulated nonresonant background with variable normalization plus a Gaussian with 
variable position, width, and magnitude. The width of the observed 0^ peak is always 
consistent with the experimental resolution of am — 3.5 MeV. Indeed, cutting on cos 0^® 
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Rejected cos 0g?^® interval 

mo (MeV) 

Signal (ev) 

s/Vb 

— In L 

xVndf 
(signal fit) 

— InL 

xVndf 
(null fit) 

2A In L 

Stat. 

sign. 

None 

1538 ± 1 

74.9 ± 14.5 

6.8 

31.5 

64.5/62 

46.6 

91.3/64 

30.1 

5.1a 

-0.6 < COS 05 P® < +0.2 

1538 ± 1 

64.1 ± 11.3 

8.6 

35.0 

59.6/62 

57.2 

91.9/64 

44.4 

6.3a 

-0.7 < cos0^® < +0.3 

1538 ± 1 

59.3 + 10.3 

9.3 

30.9 

54.4/62 

55.2 

87.5/64 

48.5 

6.6a 

-0.75 < cos 05 ^® < +0.35 

1538 ± 1 

55.8 + 9.8 

9.4 

29.3 

49.5/62 

54.4 

81.8/64 

50.2 

6.8a 


Table 1: The results of the fits of the pK^ mass spectra under the selections pr < 300 
MeV and px > 0, in which the Gaussian width of the signal has been constrained to the 
simulated resolution of am = 3.5 MeV. The statistical significance of the signal, estimated 
using the method of maximum likelihood, is shown in the rightmost column. Also shown 
is the “naive” estimate of the statistical significance S/y/B, where the signal S and the 
background B are derived from the signal hypothesis alone over the 90% area of the 
Gaussian. 

is seen to result in a dramatic increase of the signal-to-background ratio as compared to 
the pK^ mass spectra of Fig. [3l 

In order to reduce the number of free parameters, the width of the peak is con¬ 
strained to the simulated value of am = 3.5 MeV when estimating the statistical signif¬ 
icance of the signal. The results of the constrained hts of the pK^ mass spectra under 
the selections pt < 300 MeV and > 0 are shown in Table [1] Also shown for each 
fit is the difference between the log-likelihood values for the signal and null hypotheses, 
—2A In L. For the constrained fits, the numbers of degrees of freedom for the signal and 
null hypotheses differ by Andf = 2. The statistical significance of the signal is estimated 
using the value of degree of freedom which corresponds to the same p-value 

as = —2A In L for two degrees of freedom. Rejecting the central values of cos 0^® 
is seen to boost the statistical signihcance of the 0+ signal from 5.la up to 6.8a. The 
“naive” estimate of the statistical signihcance reaches S/yfB = 9.4a, where the signal 
S and the background B have been derived from the signal ht alone over the 90% area 
of the Gaussian. That the signihcance of the signal is substantially increased by an 
asymmetric cos 0™® cut a posteriori indicates that both the quadratic and linear terms 
contribute to the cos 0^® distribution. 
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Figure 7: Under the selections pt < 300 MeV and pi > 0, beam momentum p[K^) 
plotted versus the pK^ effective mass (a). The effect of rejecting the events with —0.75 < 
cos0^® < +0.35 is shown in (b). Either scatter plot is htted to the corresponding 
simulated distribution with variable normalization plus a two-dimensional Gaussian. The 
width of the Gaussian in m{pK^) and the correlation parameter p have been constrained 
to the experimental mass resolution am = 3.5 MeV and to zero, respectively. The ellipses 
show the 90% areas of the Gaussians. 


And finally, beam momentum p{K~^) is plotted versus the pK^ effective mass inder 
the selections pt < 300 MeV and > 0 in Fig. Ha), and upon rejecting the events with 
—0.75 < cos0^® < +0.35 — in Fig. Hb)- The latter selection results in a distinct 0+ 
signal at expected values of m{pK^) and p{K^). Either scatter plot is then htted to the 
corresponding simulated distribution with variable normalization plus a two-dimensional 
Gaussian. The width of the Gaussian in m{pK^) is constrained to the experimental mass 
resolution of am = 3.5 MeV, and the correlation parameter p is constrained to zero as 
physically expected for formation of a narrow resonance (as the observed mass should not 
depend on beam momentum). Either scatter plot has also been htted to the background 
form alone (not shown). For the hts in Figs. H^) andH^), we have —2AlnL = 49.7 
and 61.9 for Andf = 4. The statistical signihcance of the signal is again estimated 
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using the value of for one degree of freedom which corresponds to the same p-value 
as = —2AlnL for four degrees of freedom. Thereby, we obtain that cutting away the 
COS0™® region of —0.75 < cos0^® < +0.35 allows to increase the statistical signihcance 
of the 0’'' —?• pK^ signal from 6.2cr up to 7.la. 

In summary, using the data on the charge-exchange reaction K^n —)■ pK^ on a 
bound neutron, we have analyzed the dependence of the background-subtracted 0’'' —)■ 
pK^ signal on the emission angle in the pK^ rest frame, 0^7 In order to describe the 
observed cos 0™® distribution, invoking the interference between the nonresonant s-wave 
and the 0’''-mediated higher-wave contributions to the amplitude of the charge-exchange 
reaction is required at a level of 2.8a. The spin-parity assignment of 1/2“ for the 0’'' 
baryon is ruled out at a statistical level of 2.9a. A physically-meaningful selection in 
cos 0™® based on the observed angular dependence of the background-subtracted 0’'' —)■ 
pK^ signal allows to boost the statistical signihcance of the signal up to 6.8a for the one¬ 
dimensional pK^ mass spectrum, and 7.1a for the scatter plot in m{pK^) and p{K~^). This 
is far in excess of previously reported signals P[I21[I7] and renders the 0’'' existence more 
credible [21]. A high-statistics investigation of the charge-exchange reaction K^n —)■ pK^ 
is needed for hnally proving the 0''' existence and hxing its quantum numbers. 

Instructive discussions with professor Ya. Azimov of the St. Petersburg Institute 
of Nuclear Physics are gratefully acknowledged. 
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